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The  effects  of  self-generated  magnetic  field  in  laser 
produced  plasmas  on  the  parametric  decay  of  an  extraordinary 
electromagnetic  wave  into  two  upper  hybrid  plasmons  is 
examined  for  arbitrary  magnetic  field  intensity  and  arbitrary 

ratio  k/kn.  Due  to  the  presence  of  magnetic  field,  the  linear 
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Landau  damping  is  greatly  reduced  and  the  spectrum  of  unstable 
modes  is  significantly  modified  for  kX^>0.2. 
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I . INTRODUCTION 


A dc  magnetic  field  of  several  megaGauss  is  generated  spontaneously 
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near  the  critical  density  layer  of  a laser  irradiated  plasma. 

Because  of  the  plasma  expansion,  the  self-generated  magnetic  field 

spreads  out  to  the  underdense  regions  of  the  plasma.  It  has  been 
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measured  experimentally  and  in  computer  simulations.  In  the  present 
paper  we  study  the  effects  of  this  magnetic  field  on  the  parametric 
decay  of  the  laser  light  into  two  plasma  waves  near  the  quarter  critical 
density. 

The  nonlinear  process  of  an  electromagnetic  wave  decaying  into 

two  plasma  waves  in  a homogeneous  and  unmagnetized  plasma  was  first 
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studied  by  Goldman  and  Jackson.  Goldman  analyzed  the  process  by 

examining  the  Green  function  for  Poisson's  equation  when  the  pump 

induced  energy  of  the  particles  is  small  compared  to  the  thermal 

energy.  He  showed  the  existence  of  the  instability  by  keeping  the 

pump  wavenumber  finite  since  the  dominant  part  of  the  nonlinear 

susceptibility  is  proportional  to  it,  and  found  the  threshold  for 

the  instability.  Jackson  used  linearized  Vlasov  equation  allowing 

the  pump  Intensity  to  be  well  above  threshold.  He  showed  that  the 

system  is  stable  in  the  dipole  approximation;  however,  if  the  pump 

wavenumber  is  finite  then  the  most  unstable  perturbations  are  those 

for  which  the  wave  vector  of  the  decay  wave  lies  in  the  plane 

determined  by  the  propagation  and  polarization  vectors  of  the  pump 

wave  and  bisects  the  right  angle  between  them.  He  found  a threshold 

condition  one  order  of  magnitude  higher  than  that  calculated  by 

Goldman.  More  recently,  there  has  been  renewed  interest  in  the  parametric 
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decay  of  laser  radiation  into  two  plasmons  in  an  inhomogeneous  plasma. 
Rosenbluth*1  used  WKB  approximation  to  derive  the  threshold  condition 
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for  the  growth  of  the  decay  waves  in  order  to  overcome  the  convective 

12 

loss  out  of  the  three-wave  resonance  region.  Lee  and  Kaw  showed 

13 

the  absolute  nature  of  this  parametric  instability.  Liu  and  Rosenbluth 

used  an  alternative  method  to  analyze  the  instability  instead  of  the 

usual  WKB  approximation  and  found  the  growth  rate,  the  absolute  instability 

condition,  the  threshold  condition  imposed  by  plasma  inhomogeneity, 

and  the  saturation  level  of  the  plasma  waves  due  to  pump  depletion. 

14 

Schuss  considered  the  problem  of  an  electromagnetic  wave  obliquely 
incident  on  the  density  gradient  and  found  that  the  threshold  of  the 
absolute  instability  decreases  as  the  angle  between  the  density  gradient 
and  the  propagation  vector  approaches  90°.  Experimentally,  the  decay  of 
electromagnetic  wave  into  two  plasmons  has  been  observed  through 
measurements  of  plasma  emissions  at  the  three-halves  harmonic 

of  the  incident  laser  frequency  from  the  quarter  critical  density 
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layer. 

In  this  work  we  investigate  further  the  linear  instability 
properties  by  including  the  effects  of  the  self-generated  dc  magnetic 
field  on  the  decay  of  laser  light  into  two  upper  hybrid  plasmons. 

These  effects  might  be  relevant  since  the  magnetic  field  changes 
the  plasma  dispersive  properties,  for  instance,  if  the  wave  vectors 
of  the  decay  waves  are  in  the  plane  perpendicular  to  the  magnetic  field 
the  growth  rate  is  finite  even  for  kAD=l,  is  the  Debye  length, 
contrary  to  the  unmagnetized  plasma  where  the  decay  waves  are  heavily 
Landau  damped  for  kXjj>0.2  and  the  instability  is  turned  off. 

In  Sec.  II,  we  derive  the  system  of  nonlinear  coupled  equations 
describing  the  decay  of  an  extraordinary  mode  into  two  upper  hybrid 
waves.  The  dispersion  relation  is  solved  to  find  the  growth  rate, 
in  Sec.  III.  In  Sec.  IV,  we  present  the  conclusions  and  the  numerical 


’J  jl 


4 


results  on  the  variation  of  growth  rate  with  magnetic  field  intensity, 
pump  and  decay  wave  wavenumbers,  and  the  angle  between  these  wave  vectors. 


II.  NONLINEAR  COUPLED  EQUATIONS 


Consider  an  electromagnetic  wave 


EQ (x , t ) = -y  [exp(i^o*x-i(D0t)+c.c.  ] 


in  a magnetized  plasma  with  the  uniform  magnetic  field  B^  along  the 
z-direction.  Assuming  the  electromagnetic  pump  to  be  an  extraordinary 
mode,  we  find  electron  oscillations  along  and  perpendicular  to  the 
direction  of  propagation  with  velocities  to  be 
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where  « =eB„/mc,  n^  is  the  particle  density  of  the  unperturbed  electron 
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fluid,  e is  the  electron  charge,  and  c is  the  speed  of  light.  The 
electrostatic  component  of  the  electric  field  can  be  written  in  terms 
of  the  electromagnetic  component  as 
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The  pump  wave  magnetic  field  la  given  by 


In  - {.«*„ 
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The  perturbed  density  fluctuations  for  the  decay  waves  are 
obtained  from  the  equation  of  continuity  which,  after  Fourier 
analysis,  gives 
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where  tt  and  to  “uj-(Oq  according  to  the  resonance  conditions. 


The  anti-Stokes  component  is  considered  off-resonant  for  this 
parametric  process.  The  perturbed  velocities  v and  v are  calculated 
from  the  equations  of  motion 
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where  T is  the  electron  thermal  energy,  and  ? and  are  the  perturbed 
electric  fields.  Fourier  analyzing  Eqs.  (11)  and  (12)  for  (£,u>)  and 
(It  ,o)  ),  respectively,  we  get  the  equations  for  the  components  of  v and 
v which  together  with  Eqs.  (2) -(4)  are  substituted  into  Eqs.  (9)  and 
(10)  to  get  the  expressions  for  the  perturbed  density  oscillations 
for  the  two  decay  waves.  Substituting  these  expressions  for  n and  n_ 
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into  Poisson' 8 equations  for  the  perturbed  electrostatic  potentials 
4 and  4 we  get  the  system  of  coupled  equations 
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where  n and  n are  the  nonlinear  contributions  to  the  density 
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and  M is  the  ion  mass.  The  wavenumbers  k„  and  k_|(  , kA  and  k_A 
refer  to  the  components  of  £ and  parallel  and  perpendicular 
to  the  magnetic  field,  respectively.  The  ions  are  considered 
cold  and  unmagnetized  since  the  frequencies  of  the  pump  and  decay 
waves  are  much  larger  than  the  ion  cyclotron  frequency.  Equations  (13) 
and  (14)  can  be  rewritten  as 
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where 
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and  z*z/z.  and  a"!  are  obtained  from  and  cij,  , respectively,  by 
making  the  following  interchanges 

ai“«I I (£,u>)-*-*-(ic_ ,u»_)  , (£0,w0)-~.(-£0,-a)0) , 
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and 

a",  -a],  [(£,w)-*-*-Oi_,w_)  , (kQ.w^-'-K-kg.-ujg) , vQ->v*]  . (23) 

Equations  (18)  and  (19)  comprise  the  system  of  coupled  equations 
describing  the  parametric  decay  of  an  extraordinary  electromagnetic 
wave  into  two  electrostatic  waves.  The  general  dispersion  relation 
is  obtained  from  them,  straightforwardly, 

ee_=(a]+a[(  )(a^+a']'|  ) . (24) 

Equation  (24)  allows  us  to  study  the  decay  of  an  extraordinary 
electromagnetic  wave  into 

(i)  two  upper  hybrid  waves, 

(ii)  an  upper  hybrid  and  a lower  hybrid  wave,  or 

(iii)  two  lower  hybrid  waves. 

The  decay  (i)  occurs  at  the  quarter  critical  density  while  channel  (ii) 
at  the  critical  density  layer.  In  the  present  paper  we  restrict  ourselves 
to  channel  (i). 


III.  GROWTH  RATE 

For  the  decay  of  an  electromagnetic  wave  into  two  upper  hybrid 
plasmons,  the  linear  dispersion  relations  (15)  and  (16),  in  the  limit 
k()<<  k,  become 
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If  we  set  ft  -0  in  Eqs.  (27)  and  (28)  we  get  the  same  equations  as  Liu 
e 

and  Rosenbluth  in  their  limit  L-*»  where  L is  the  density  scale  length. 
For  laser  fusion  parameters,  fte/tOp<l,  Eqs.  (27)  and  (28)  give  the 
dispersion  relation 

[^-(uV+kV)  i [(A((oW+kV)  ] 
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The  growth  rate  is  found,  from  Eq.  (29),  to  be 
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IV.  NUMERICAL  RESULTS  AND  CONCLUSIONS 

The  growth  rate  decreases  slowly  with  the  Increasing  magnetic  field 
intensity,  as  seen  from  Fig.  1.  The  maximum  growth  rate  is  Ymax  * 

k0|v0|  which  holds  for  fie=0,  result  that  agrees  with  Liu  and  Rosenbluth. 
The  linear  Landau  damping  for  the  upper  hybrid  waves  with  finite  k)( 
is  given  by 
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and  p£  is  the  electron  Larmor  radius.  From  Eq.  (31)  we  infer  that  if 

k and  k are  in  the  plane  perpendicular  to  the  magnetic  field  the  upper 

hybrid  decay  waves  cannot  resonate  with  the  electrons  and,  therefore, 

the  linear  Landau  damping  rate  vanishes  in  this  case.  It  means  that 

the  growth  rate  spectrum  is  significantly  modified  for  kXp>0.2  as 

compared  to  the  unmagnetized  plasma  where  we  would  expect  the  decay 

waves  to  be  strongly  Landau  damped.  Figure  2 shows  the  growth 

rate  as  a function  of  kA^  for  various  pump  powers.  The  growth  rate 

is  finite  for  all  values  of  kA^  due  to  the  absence  of  collisionless 

damping,  except  for  k=kg/2cos<(>  for  cos<f>>0  when  it  vanishes,  d>  is  the 

angle  between  k^  and  ic.  The  growth  rate  is  always  positive  for  cos<jx0. 

According  to  Fig.  1,  for  laser  fusion  parameters,  ft  « w , the  growth 

e p 

rates  for  magnetized  and  unmagnetized  plasmas  have  the  same  magnitude 
for  kAjj<0.2.  However,  for  kAD>0.2  they  are  substantially  different 
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due  to  the  inclusion  of  Landau  damping,  as  can  be  seen  from  Fig.  3. 
For  instance,  for  laser  powers  up  to  lO^^W/cm^  the  growth  rate 
vanishes  for  kXD>0.23  for  an  unmagnetized  plasma  process  but  it  is 
finite  if  the  magnetic  field  is  present.  It  is  also  possible  to 
see  from  Fig.  2 that  as  cos<(>  increases  the  values  of  kX^  decrease  for 
the  maximum  growth  rate.  This  is  better  seen  from  Fig.  4 where  <p  is 
plotted  versus  kX^  for  a constant  growth  rate.  Figure  5 exhibits  the 
variation  of  the  growth  rate  with  cost)).  The  growth  rate  vanishes  for 
cos<fi=ko/2k  if  kg<2k.  For  kg>2k  the  growth  rate  is  finite  for  all 
values  of  cos<j>.  Figure  6 shows  the  proportionality  of  the  growth 
rate  with  k^X^. 
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FIGURE  CAPTIONS 

Fig.  1 Variation  of  the  growth  rate  Y*d/|vq|  with  the  electron  cyclotron 

frequency  for  kgAjj=0.05  (0.32  keV) , kAp=0.002,  and  cos$=0.0. 

Fig.  2 Growth  rate  Y*d/|vqI  expressed  as  a function  of  kAp  for  kgAp-0.05, 

n2/o)2=0.01,  cos<J>=0.0  ( ),  0.6  ( ),  -0.6  (-0-0—),  0.9  (———), 

e p 

and  -0.9  (-0-0-). 

Fig.  3 Variation  of  the  growth  rate  y/Wp  with  kip  for  kgAD=0.1  (1.275  keV) , 
2 2 

0^/o)p=0.01,  for  the  following  pump  powers:  (i)  |vg|/ve=1.17 

(1015W/cm2),  cos<|>=0.6  ( ),  -0.6  (-KM)-) ; (ii)  |vQ|/ve= 

3.80  (lO^W/cm2) , cos<(>=0.6  ( ),  -0.6  (■♦♦■).  The  growth 

rates  follow  curves  (-•-•-)  for  unmagnetized  plasma. 

Fig.  4 Variation  of  kAp  with  (angle  between  kg  and  1c)  for  a constant 

2 2 

growth  rate  y!d/|vq|  where  kgAp-0.05  and  Ge/u)p=0.01. 

Fig.  5 The  growth  rate  yA^/ / Vq J expressed  as  a function  of  cos$  for 

knA  =0.05,  fi2/u)2=0.0l  and  kA  =0.03  (-0-0-) , 0.15  ( ). 

u u e p u 

2 2 

Fig.  6 Variation  of  Y*d/|vq|  with  kgA^  for  =0.01  in  the  following 

cases:  (i)  kAD=0.03,  cos<t>=0.0  ( ),  0.7  (— 0-0— ) , -0.7  ( ) ; 

(ii)  kAD=0.15,  cos<f>=0.0  (-.-•-) . 


i 


C\J  o 00  CD  C\J  03 

c0lx|°A|/°\/ 

Fig.  2 


o 

<7> 


IO 


O IO 

<0 


£0lxV>C 


O 

ro 


Fig.  3 


COS  <p 


